Quantum State Tomography via Linear Regression Estimation 
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A simple yet efficient method of linear regression estimation (LRE) is presented for quantum state 
tomography. In this method, quantum state reconstruction is converted into a parameter estimation 
problem of a linear regression model and the least-squares method is employed to estimate the un- 
known parameters. The asymptotic mean squared error (MSE) bound of the estimate can be given 
analytically, which can guide one to choose optimal measurement sets. The LRE is asymptotically 
optimal in the sense that the MSE may achieve the Cramer-Rao bound asymptotically. The com- 
putational complexity of LRE is 0(d'*) where d is the dimension of the quantum state. Numerical 
examples show that LRE is much faster than maximum-likelihood estimation for quantum state 
tomography. 

PACS numbers: 03.65.Wj, 02.50.-r, 03.67.-a 



One of the essential tasks in quantum information 

Socessing is to verify the integrity of a quantum state 
. Quantum state tomography has become a standard 
technology for inferring the state of a quantum system 
through appropriate measurements and estimation [2- 
|8|. To reconstruct a quantum state, one may first per- 
form measurements on a collection of identically pre- 
pared copies of a quantum system (data collection) and 
then infer the quantum state from these measurement 
outcomes using appropriate estimation algorithms (data 
analysis). Measurement on a quantum system generally 
gives a probabilistic result and an individual measure- 
ment outcome only provides limited information on the 
state of the system, even when an ideal measurement de- 
vice is used. In principle, an infinite number of mea- 
surements are required to determine a quantum state 
precisely. However, practical quantum state tomogra- 
phy consists of only finite measurements and appropri- 
ate estimation algorithms. Hence, the choice of optimal 
measurement sets and the design of efficient estimation 
algorithms are two critical issues in quantum state to- 
mography. 

Many results have been presented for choosing optimal 
measurement sets to increase the estimation accuracy 
and efficiency in quantum state tomography |9l-llll| . Sev- 
eral sound choices that can provide excellent performance 
for tomography are, for instance, tetrahedron measure- 
ment bases, cube measurement sets, and mutually un- 
biased bases [Uj. However, for most existing results, 
the optimality of a given measurement set is only veri- 
fied through numerical results [ll| . There are few meth- 



ods that can explicitly give an estimation error bound 



12l |13| , which is essential to evaluate the optimality of 
a measurement set [lj-[l6| and the appropriateness of an 
estimation method. 

For estimation algorithms, several useful methods in- 
cluding maximum-likelihood estimation (MLE) [a llTl - 
^, Bayesian mean estimation (BME) [1 [2l|, Q and 
least-squares (LS) inversion |23j have been proposed for 
quantum state reconstruction. The MLE method simply 
chooses the state estimate that gives the observed results 
with the highest probability. This method is asymptoti- 
cally optimal in the sense that the estimation error can 
asymptotically achieve the Cramer- Rao bound. However, 
MLE usually involves solving a large number of nonlinear 
equations where their solutions are notoriously difficult 
to obtain and often not unique. Recently, an efficient 
method has been proposed for computing the maximum- 
likelihood quantum state from measurements with ad- 
ditive Gaussian noise, but this method is not general 
[20j . Compared to MLE, BME can always give a unique 
state estimate, since it constructs a state from an in- 
tegral averaging over all possible quantum states with 
proper weights. The high computational complexity of 
this method significantly limits its application. The LS 
inversion method can be applied when measurable quan- 
tities exist that are linearly related to all density ma- 
trix elements of the quantum state being reconstructed 
[23|. However, the estimation result may be a nonphys- 
ical state and the mean squared error (MSE) bound of 
the estimate cannot be determined analytically. 

In this Letter, we present a new linear regression es- 



timation (LRE) method for quantum state tomography 
that can identify optimal measurement sets and recon- 
struct a quantum state efficiently. We first convert the 
quantum state reconstruction into a parameter estima- 
tion problem of a linear regression model 2J]. Next, we 
employ an LS algorithm to estimate the unknown pa- 
rameters. By exploring the statistic property of the lin- 
ear regression model, we are able to determine the MSE 
bound and demonstrate its asymptotic optimality in the 
minimal estimation error sense. The positivity of the 
reconstructed state can be guaranteed by an additional 
least-squares minimization problem. The total compu- 
tational complexity is 0{d'^) where d is the dimension 
of the quantum state. The efficiency of the method is 
demonstrated by examples on qubit systems. 

Linear regression model. We first convert the quantum 
state tomography problem into a parameter estimation 
problem of a linear regression model. Suppose the di- 
mension of the Hilbert space Ti of the system of interest 

i2_-| 

is d, and {^i}i^Q is a complete basis set of orthonormal 
operators on the corresponding Liouville space, namely, 



Tr(^I^ 
joint and u^j 



Jij, 



where f denotes the Hermitian ad- 
6ij is the Kronecker function. Without loss 
of generality, let fli = Q] and Qq ~ {l/d)^I, such that 
the other bases are traceless. That is Tr(f2i) = 0, for 
1 = 1, 2, • • • , d^ — 1. The quantum state p to be recon- 
structed may be parameterized as 



/ 



rf^-i 



P = 



J2 ©^"^ 



(1) 



where Qi = Tr{pi^i). Given a set of measurement bases 
{|v[/)(\]/|(")}^£^^ each |*)(vl/|(") can be parameterized un- 

der the bases {^i}^^'^ as 



d^-i 



i*)(*i^"^ = ^+E^i"^^- 



where V,|"^ = Tr(|1')(*|(")r2,). 

When one performs measurements with measurement 
set {|^)(*|'"-'}^=i on a collection of identically prepared 
copies of a quantum system (with state p), the probabil- 
ity to obtain the result of \^'){^\'^"-'> is 



Tr(|*)(*|("V) 



d ^-^ d 

i=\ 

(2) 
Assume that the total number of experiments is iV and 
N /M experiments are performed on N/M identically pre- 
pared copies of a quantum system for each measurement 
basis |\I>)(^|'"'. Denote the corresponding outcomes as 
Xi , • ■ • ,x^,j^j, which are independent and identically 



i-i , 



distributed. Let Pn 



N/M 



"n/m 



and e„ = p„ -p„. 



According to the central limit theorem (25| , e„ converges 
in distribution to a normal distribution with mean and 

_ 2 

variance ^^/^.f ■ Using ([2]), we have the linear regression 
equations 



eTxJ/(n) 



(3) 



Note that the variance of e„ is asymptotically ^Zr/fp ■ 
^i Pn = Ij we have already reconstructed the state as 
v[/^(\I/|("); if Pjj = 0, we should choose the following 
measurement basis from the orthogonal complementary 
space of |\1')(^|'^"^ Pn, d and ^'■"■' are all available for 
n = 1, • • • , M , while e„ may be considered as the ob- 
servation noise. Hence, the problem of quantum state 
tomography is converted into the estimation problem of 
the unknown vector via these linear regression equa- 
tions. 

Now, we define the MSE as ETr(/5 — p)^ where p is an 
estimate of the quantum state p based on the measure- 
ment outcomes and E(-) denotes the expectation on all 
possible measurement outcomes. For a fixed tomogra- 
phy method, ETr(/) — p)^ depends on the state p to be 
reconstructed and the chosen measurement bases. From 
a practical viewpoint, the optimality of a chosen set of 
measurement bases may rely upon a prior information 
but should not depend on the unknown quantum state 
to be reconstructed. In this Letter, no a prior assump- 
tion is made on the state p to be reconstructed. Given 
a fixed tomography method, we use the maximum MSE 
bound for all possible states (i.e., suppETr(p — p)^) as 
the index to evaluate the performance of a chosen set of 
measurement bases. Hence, it is necessary to consider 
the worst case by enlarging the variance of the obser- 
vation noise e„ in each linear regression equation. As 
a consequence, {e„}*£]^ may be treated as a set of inde- 
pendent identically distributed variables with asymptotic 
normal distribution N(0, j^)- Another advantage of this 
treatment is that the effect of some other noises can be 
absorbed in the enlarged variance. 

Asymptotic properties of the LS estimate. To give an 
estimate with high accuracy and low computational com- 
plexity, we employ the LS method, where the basic idea 
is to find an estimate Ols such that 



Qls 



M 

argmin \ {p„ 



d ' ' 



where Q is an estimate of Q. The LS solution is 

M 



eLs = pE*^"^(p„-^), 



n=l 



where 



M 



P = (y^^(")^(«) )-i^ 



(4) 



(5) 



If the measurement bases {|^)(^|''"''}*Li a-rc informa- 
tionally com.plete or overcomplcte such that P is invert- 
ible, then as A'' goes to infinity, the estimate Qls has the 
following properties |24| : 

1. Qls is asymptotically unbiased; 

2. Qls is asymptotically a maximum- likelihood es- 
timate, and the estimation error can asymptotically 
achieve the Cramer-Rao bound; 

3. The z-th element of Qls is asymptotically the min- 
imum variance unbiased estimation of Qi ; 

4. The MSE of Qls H^ls - QV{Qls - Q) is asymp- 
totically ^TrP, i.e., iiTr(Ef=i*(")*(")^)"^ 

Positivity and computational complexity. Based on the 
solution Qls obtained from ^, we can obtain a Her- 
mitian matrix fi with Tr/i ~ 1 using ([T]). However, /i 
may have negative eigenvalues and be nonphysical due 
to the randomness of measurement results. In this sense, 
jl is called pscudo linear regression estimation (PLRE) of 
state p. A good method of pulling jl back to a physical 
state can reduce the MSE. In this Letter, the physical 
estimate p is chosen to be the closest density matrix to p 
under the matrix 2-norm. In standard state reconstruc- 
tion algorithms, this task is computationally intensive 
[20|. However, we can employ the fast algorithm in [20[ 
with computational complexity 0{d^) to solve this prob- 
lem since we have obtained a Hcrmitian estimate p with 
Tip = 1. 

Since an informationally complete measurement set 
{|*)(vl/|(")}^^£^ requires M being O(d^), the computa- 
tional complexity of ^ and ^ is 0{d'^). Although the 
computational complexity of ^ is generally 0{d^), ^ 
can be computed off-line before the experiment once the 
measurement set is determined. Hence, the total compu- 
tational complexity of LRE after the data have been col- 
lected is 0{d'^). It is worth pointing out that for n-qubit 

systems, X]n=i ^^"'^VE'^"'^ is diagonal for many preferred 
measurement sets such as tetrahedron and cube measure- 
ment sets. Fig. 1 compares the run time of our algorithm 
with that of a traditional MLE algorithm. Since the max- 
imum MSE could reach 2 for the worst estimate, it is clear 
that our algorithm LRE is much more efficient than MLE 
with a small amount of accuracy sacrificed. 

Optimality of measurement bases. One of the advan- 
tages of LRE is that the MSE bound can be given as 
3^Tr(X;*ii 1'(")*(")^)-i, which is dependant explicitly 
upon the measurement bases. Recall that an optimal 
set of measurement bases should minimize sup ETr(p — 

p)2 (i.e., i^Tr(X;*ii*("^*^"^^)"^), under all possible 
states. In principle, an optimal set {|^)(^|'"''}^Li can 
be chosen by solving the following optimization problem: 



Minimize Tr(^„ii ^(")^(") ) 



s.t. *(") *("' 



d-l 

d ' 



for n — 1, 



M. 




FIG. 1: (color online) The run time and MSE of LRE and 
MLE for n-qubit state i|GHZ>(GHZ| + ^^, where |GHZ> = 
-^i\H^") + \V^")). The realization of MLE used the iterative 
method in Q]. The measurement bases are from the n-qubit 
cube measurement set and the resource is A'' = 3^ x 4". LRE 
is much more efficient than MLE witli a small amount of 
accuracy sacrificed since the maximum MSE could reach 2 for 
the worst estimate. All timings were performed in MATLAB 
on the computer with 4 cores of 3GHz Intel 15-2320 CPUs. 



way by employing appropriate algorithms though it may 
be computationally intensive [26[. 

With the help of the analytical MSE bound, we can 
ascertain which one is optimal among the available mea- 
surement sets. This is shown when we prove the opti- 
mality of several typical sets of measurement bases for 
2-qubit systems below. 

For 2-qubit systems, it is convenient to chose fli = 
-y=ai ® -i=cr„j, where i = AI + m\l^ m = 0, I, 2, 3; (Jq — 



^2x2i CTl — 



1 \ _ ( -i 

1 o)^"'^ [i 

The MSE bound of 2-qubit states is 



CT3 = 



1 
-I 



M 

4iV 



M 



TrP 



— Tr(^*("%(")^^ 

n=l 



47V 



Now we minimize this MSE bound or equivalently mini- 
mize Tr(P). Denote the eigenvalues of P~^ as Ai > A2 > 



> A15. Since we have V'o 



E!=0 V^: 



15 



i. the 



subproblem is converted into minimizing X)i=i T"' sub- 
ject to J2i=i '^i — J^^- It '^''^^ be proven that J2i=i X" 



4 
reaches its minimum ^^ when Ai 



5 ~ 20- 
^^TrP is ^. 



The optimization problem can be solved in an off-line 



300 „,i,„„ \ . — . . . ^ x-^^ ^ K 

Hence, the minimum of the MSE bound aj^ -•-•^^ ^^^ « 
This MSE bound can be reached by using the mutually 
unbiased measurement bases [9|. 

If only local measurements can be performed, i.e.. 



"S _s 
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FIG. 2: (color online) Mean squared error (MSB) for Werner 
states [231 with q (varying from to 1) and different numbers 
of copies A^. The cube measurement set is used, where the 
MSE bound is ^. It can be seen that the MSE of PLRE is 
almost unchanged for q £ [0, 1], and is larger than the MSE 
of LRE. For states which are not highly pure, the MSE of 
LRE is almost the same as that of MLE. As the purity of the 
states increases, the MSE of MLE becomes smaller than that 
of LRE. 



|$)(^|(") = |^)(*|(".i) |*)(^|(".2)^ n = 1, ••• , M, 
where |*)(\l/|("'i) and |*)(\?|("'2) can be parameterized 



as |*)(^| 



(n,fc) _ 



E/=0 ^' 



t. 



— V'; X V'm ' , where i 



.,(«,fc)_^ fc ^ 1^ 2. And we have 
Al+m. Due to additional 



V2' 



constraints tAq ' 71' ^i=o^i 



(n,k)' 



1, 



1, for fc = 1, 2, 
M, the subproblem of minimizing the 

15 



MSE bound can be converted into minimizing X]i=i T"; 



(iii) E 



15 



subject to (i) ^ti ^» > 3^^; (ii) ELi ^^ > W- 
It can be proven that Ei=i T" 
reaches its minimum -^ when Ai = • • ■ = Ag = |^, 



lA, = |Af. 



At 



A 



15 



41 . Hence, the minimum of the MSE 

36 ' 



bound 3^TrP is ^. This MSE bound can be reached by 
using the 2-qubit cube or tetrahedron measurement set 

i- 

Fig. 2 shows the dep endant relationships of the MSEs 
for Werner states 27| on q (varying from to 1) and 
different number of copies N using the cube measurement 
bases 0. The fact that the MSE of PLRE is larger than 
that of LRE demonstrates that the process of pulling /t 
back to a physical state further reduces the estimation 
error. Note that there are small gaps between the MSE 
bound 22 and the MSE of PLRE or LRE. This is because 
ideal cases are considered in our simulation. If other noise 
is taken into consideration, the MSEs of PLRE and LRE 



could touch the bound ^. 

Discussions and conclusions. In the LRE method, data 
collection is achieved by performing measurements on 
quantum systems with given measurement bases. This 
process can also be accomplished by considering the 
evolution of quantum systems with fewer measurement 
bases. For example, suppose only one observable a is 
given, and the system evolves according to a unitary 
group {Ut\- At a given time i, 

{at) = TT{UHt)aU{t)p) = Tr^p). 

Suppose one measures the observable a at time t {t ^ 
1, • ■ • , M) on m identically prepared copies of a 
quantum system. Denote the obtained outcomes as 
• • , a^^, and their algebraic average as at — 

""^ . Note that ct* , • • • , a\^^ are independent and 
identically distributed. According to the central limit 
theorem [25|, Ct^ at — {at) converges in distribution to a 
normal distribution with mean and variance ■^— ^-^ — . 

m. 

We have the following linear regression equations 






CTt = Tr(o-tp) + et, t = 1, 



M, 



which are similar to ([3|). Hence, we can use the proposed 
LRE method to accomplish quantum state tomography. 

The LRE method can also be extended to reconstruct 
quantum states with a prior information 12|, |28|, |29| or 
states of open quantum systems. Actually, LRE can be 
applied whenever there are measurable quantities that 
are linearly related to all density matrix elements of the 
quantum system under consideration. 

In conclusion, an efficient method of linear regression 
estimation has been presented for quantum state tomog- 
raphy. We have analytically provided the MSE bound 
of the estimation, which explicitly depends on the mea- 
surement bases and can assist to identify optimal mea- 
surement sets. The computational complexity of LRE is 
0[d^), which is much lower than that of MLE and BME. 
It is possible to extend other results from classical linear 
model theory to estimation problems of quantum states. 
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